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<0 Abstract 



We propose a phenomenological model that describes counterflow and drag experiments with quantum Hall bilayers in a 
vt = 1 state. We consider the system consisting of statistically distributed areas with local total filling factors vti > 1 
5-h and vti < 1- The excess or deficit of electrons in a given area results in an appearance of vortex excitations. The vortices 
in quantum Hall bilayers are charged. They are responsible for a decay of the exciton supercurrent, and, at the same 
time, contribute to the conductivity directly. The experimental temperature dependence of the counterflow and drive 
resistivities is described under accounting viscous forces applied to vortices that are the exponentially increase functions 
of the inverse temperature. The presence of defect areas where the interlayer phase coherence is destroyed completely 
pan result in an essential negative longitudinal drag resistivity as well as in a counterflow Hall resistivity. 
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1. Introduction 



Theoretical prediction on superfluidity of bound 
electron-hole pairs in electron-hole bilayers [IrQ and, es- 
pecially, in quantum Hall bilayers [5l4l3j has stimulated 
experimental study of the counterflow transport in such 
systems. In these experiments [14H19I ] electrical current is 
passed through one layer in a given direction and returned 
to the source through the other layer in the opposite direc- 
tion. The counterflow currents in the adjacent layers can 
be provided by excitons that consist of an electron belong- 
ing to one layer and a hole belonging to the other layer. 
One can expect that at low temperature the gas of such 
excitons becomes superfluid and the counterflow current 
may flow without dissipation. 

In the counterflow experiments 14 - 19| bilayer quantum 
Hall systems with the total filling factor vt = 1 are used. 
At such a filling the number of electrons in one layer coin- 
cides with the number of empty quantum states (holes) in 
the lowest Landau level in the other layer. At rather small 
interlayer distances [2(| the ground state of such a sys- 
tem is the BCS-like state with electron-hole pairing. The 
pairing is caused by the Coulomb attraction. The BCS- 
like state can also be considered as a superfluid state of 
interlayer excitons. 

In spite of expectation to reach at low temperature zero 
electrical resistance, counterflow experiments in quantum 
Hall bilayers give some other results. At all temperatures 
a finite longitudinal resistance is registered. It decreases 
exponentially under lowering of temperature, but it does 
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not vanish completely. Moreover, the counterflow longitu- 
dinal resistivity p£,f is even higher than the longitudinal 
resistivity measured in the parallel current geometry p\ x 
(in which the electrical currents in the adjacent layers are 
equal in value and flow in the same direction). 

The excitons cannot give a contribution into the con- 
ductivity in the parallel current geometry. Therefore, the 
parallel current conductivity o~\ x should be much smaller 
than the counterflow conductivity a xx . The relation be- 
tween the longitudinal resistivities is another one. At zero 
or small Hall conductivity the longitudinal resistivity is 
in inverse proportion with the longitudinal conductivity. 
But for the conductivity tensor with a large Hall compo- 
nent the longitudinal resistivity is in direct proportion with 
the longitudinal conductivity. The first case is realized for 
the counterflow geometry and p xx ~ . The parallel 

current geometry corresponds to the second situation and 

pxx ~ 0~xx/{&xy) • 

Exciton superfluidity can be imperfect due to phase slips 
caused by motion of vortices across the flow. In that case 
the counterflow resistivity is nonzero and, in principle, it 
can be larger than the parallel current resistivity. The 
smallness of p\ x can be accounted, for instance, for a small 
density of carriers that contribute into cr| x . 

The specifics of the exciton superfluidity in bilayer quan- 
tum Hall systems is that the vortices carry electrical 
charges. Such vortices not only bring to a decay of su- 
perflow, but also give a direct contribution into conduc- 
tivity. As was argued in [2l(, due to such a two- fold role 
of vortices, the counterflow and parallel current longitudi- 
nal resistivities coincide with each other (in the absence of 
other essential factors that may contribute to the conduc- 
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tivity) . 



The role of vortices in transport behavior of quantum 
Hall bilayers was considered in Ref. 
ation 



22] 



The consider- 

is based on the coherence network model j23[. 
In this model superfluid excitons flow in a network formed 
by narrow links, and the dissipation is caused by vortices 
that cross these links. In [22J the draegeometry experi- 
ments are analyzed. In this geometry the current 
flows through one layer and the resistivities are measured 
in both layers. The resistivities in the active (drive) layer 
" ' ,e and in the passive (drag) layer p^™ 9 are connected 

.11 i „CF\ 



p: 

with p'^p and by the relations p' a8 



drive 



i drag 
and P a t ~ 



- 0>i(i ~ (Kp)/' 4 ' Wnere P'afi 

resistivities of the layer that is used as an active one in 
the drag geometry (in case of imbalanced bilayers p^g and 
Po$ are different for different layers). The resistivities 
Pxx V& an d Pxx demonstrate thermally activated behav- 
ior. The activation energy depends asymmetrically on the 
imbalance of electron densities of the layers 0,[l7j]- As 
was shown in [22J this behavior can be described under 
assumption that vortex mobilities are thermally activated 
quantities and their activation energies are different for 
different specie of vortices. The vortex configurations in 
strongly disordered quantum Hall bilayers were studied in 
[2~il ]. It was shown that disorder should be rather strong 
to provide vortex proliferation. 



pgf)/2, 



where 



J 



and p^p are the 



The ideas of |2Ir424j are important for the understand- 
ing of transport properties of quantum Hall bilayers. Nev- 
ertheless, the question requires further study because a 
number of essential features should be explained. First of 
all, in all experiments the drag resistivity p^™ 9 is nega- 
tive and the counterflow resistivity p^f is larger than the 
parallel current resistivity p xx . Actually, it is one of key 
points because in case of perfect exciton superfluidity p x r x a9 
should be positive and equal to p x w x lve - Then, the origin 
of the counterflow Hall resistivity remains unclear. Also, 
the absolute value of p xx observed experimentally is quite 
large and one can think about additional factors (beside 
vortices) that increase p x 
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J xx 



In this paper we present a model in which transport fea- 
tures of quantum Hall bilayers are accounted for a special 
influence of imperfectness. We assume that due to imper- 
fectness the local total filling factor deviates from unity 
and in some area Vt < l f while in other areas vt > 1< 
In the areas of the first (second) type positive-charged 
(negative-charged) vortices emerge. Besides, we imply 
that in some defect areas the interlayer phase coherence 
is destroyed completely. Using the effective medium ap- 
proach 25-2?) we compute the effective longitudinal and 
Hall resistivities for the counterflow, parallel current and 
drag geometries. The results are in good agreement with 
the experiment. 



Table 1: The vorticity a K , the electrical charge q K and the layer i K 
where the charge is concentrated, for four species of vortices 



K 


a K 


<1k 


i K 


I 


-1 


e(l-v) 


2 


2 


+1 


ev 


I 


3 


-1 


—ev 


2 


4 


+1 


-e(l-iz) 


I 



2. Vortices in exciton superfluids in quantum Hall 
bilayers 

There are four species of vortices in the system consid- 
ered [I)]. Below they are notated by the index k = 1, 2, 3, 4. 
The wave functions \k) for the states with a single vortex 
can be presented as 
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Here i is the operator of creation for an electron in the 
layer i in the state with definite angular momentum m: 
# m (r,0) = r m e-'' mf> exp(-r 2 /M 2 ), where 6 is the polar 
angle counted according to the right-hand rule with re- 
spect to the magnetic field direction, and I = \Jhcj eB is 
the magnetic length. The u and v coefficients are u = \pv 
v, where v = V\ is the filling factor for the 
v = V2 is the filling factor for the layer 2, 



and v = \/T" 
layer I, and I 

vi = 27r^ 2 n e .i, and n e ^ is the electron density in the layer 



Each vortex is characterized by distinct values of its vor- 
ticity and electrical charge. The vortex charge q K can be 
found from the corresponding wave function (TTJ) . One can 
see from (fTJ) that the charge of a given vortex is fractional 
and is concentrated in a certain layer. We notate this layer 
by i K . 

The vortex is a state with a clockwise circular electrical 
current in one layer and a counterclockwise current in the 
other layer. We use the convention that the vortex has the 
positive vorticity a K = +1 if the electrical current associ- 
ated with the vortex is a counterclockwise in the layer 1 
and clockwise in the layer 2. The vortex parameters for 
each specie are given in Table Q~] 

The vortices nucleate in pairs. The components of the 
pair have the opposite vorticities. In the absence of elec- 
tron deficit or excess the vortices in the pair have the op- 
posite electrical charges, as well. The Coulomb attrac- 
tion results in an increase of the binding energy of the 
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pair. Therefore, free vortices at vt = 1 emerge at tem- 
peratures higher than the Berezinskii-Kosterlits-Thouless 
(BKT) transition temperatur^ 

Deviation of the total filling factor from unity forces the 
appearance of charge excitations that are transformed into 
vortex pairs. In a given pair each vortex has the fractional 
charge, but the sum of the charges is integer. In areas with 
electron deficit an equal number of vortices k — 1 and 2 
is nucleated, and equal number of vortices k = 3 and 4 
emerge in areas with electron excess. In a given area the 
vortices have the same sign of charge and due to Coulomb 
repulsion the size of a bound vortex pair can be quite large. 

The energy of the pair is 



E = Ttp s In 



r e 2 v(l — v) 



sr 



(2) 



where p s is the superfluid stiffness, r is the distance be- 
tween the vortices, and a is the vortex core radius. The en- 
ergy @ is minimum at the distance = e 2 v(l — is)/irep s . 
The quantity 7 = ri,(2n v )^ 1 ^ 2 (where n v is the concentra- 
tion for one specie of vortices) yields the ratio of the vortex 
pair size to the average distance between the vortices. One 
can expect that at 7 > 1 a plasma of free vortices of op- 
posite vorticities instead of a gas of vortex pairs emerges. 
Using the mean field value of p s Q one finds that the 
condition 7=1 corresponds to the critical vortex density 
n VtC ~ 0.05no, where uq — \/2tt1 2 is the electron density 
of a completely filled Landau level. The critical density 
n„ jC is proportional to p 2 . It is known that the mean- field 
approximation overestimates the superfluid stiffness p s : it 
does not take into account quantum and thermal excita- 
tions, and the interaction with impurities that reduce the 
superfluid stiffness (see, for instance, [28l - [30j j V Therefore, 
the actual critical density can be considerable smaller. 

3. The resistivity caused by vortex motion 

Let us consider the forces that act on a vortex. The vor- 
tices are electrically charged and the electric and magnetic 
components of the Lorentz force are applied to them: 

F L = (? K E lK + q K v K x B/c, 

where is the electrical fields in the layer i and v K is the 
vortex velocity. 

Vortices carry the vorticity and in system with nonzero 
net current the Magnus force Fj.; emerges. This force can 
be obtained from the derivative of the electrical current en- 
ergy with respect to the vortex position (see, for instance, 

m 



F M = 2iTa K p s 'V(p x z. 

Here z is the unit vector directed along the magnetic field, 
and Vtp is the gradient of the phase of the order parameter 



x We do not consider the effect of vortex unbinding caused by 
electical currents, ft results in nonlinear dependence of the voltage 
on the current. 



taken far from the vortex center. The net currents in the 
layers read as j s i = — j S 2 = ep s \Jip/h. The Magnus force 
does not depend on the vortex velocity, because in the case 
considered this force is proportional to the difference of the 
net currents j s i — j S 2 . This difference is the same in the lab 
reference frame and the reference frame connected with a 
moving vortex. 

We should also take into account the viscous force 
F, ; = — r/ K v K , where rj K is the viscosity parameter. In what 
follows we will imply the vortex motion is connected with 
thermally activated hops of vortices between pinning cen- 
ters. Is this case temperature dependence of the viscosity 
parameters can be approximated as r\ K oc exp(A K /T) 

In the stationary state the resultant force applied to the 
vortex is equal to zero 

2na K ,p s X7ip X z - r] K v K + q K E lK + — v K x B = 0. (3) 

c 

Eq. (J3|) is fulfilled for each specie n with nonzero density. 

The motion of vortices across the flow results in a decay 
of the phase gradient. The rate of decay is d(V<p)/dt = 
J2 K 2irn K a K v K x z, where n K is the density of vortices of 
specie k. Antiparallel electrical fields applied to the lay- 
ers may compensate this decay. They increase the phase 
gradient at the rate d(V<p)/dt = (e//i)(Ei - E 2 ). In the 
stationary state the following condition should be satisfied 



2irn K a K v K x z + ^(Ei 
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(4) 



Having the vortex density one solves Eqs. © and (Q| 
and finds the contribution of excitons and direct contri- 
bution of vortices into conductivity. To compute the con- 
ductivity tensor one should also take into account the bare 
conductivity of Landau level. 

In the BCS-like state each electron is distributed be- 
tween two one-particle states with the same guiding center 
index and distinct layer indices. The system responds on 
the resultant electrical force Fe = — iV e (^ieEi + ^ 2 eE 2 ) 
(N e is the total number of electrons) in the same manner 
as a completely filled Landau level responds on electrical 
field. In the latter case the electron gas moves as a whole 
with the drift velocity vo = cE x B/5 2 . In case of the bi- 
layer system the expression for the drift velocity contains 
the effective field E e // = ^iEi + z^E 2 . The velocity v 
can be also found from the condition that the electrical 
and magnetic components of the Lorentz force applied to 
the whole system compensate each other, as it takes place 
for a completely filled Landau level with negligible small 
relaxation 



z^iEi + z^ 2 E 2 + -Vo x B = 0. 



(5) 



The motion of an electron gas with the velocity vo 
yields the contribution into the electrical currents <5j; = 
-en ^v . 

Thus, the electrical currents in the layer read as 
ji = -<r vi {viEi + U2E2) x z 
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+-p s Vip + q 2 n 2 v 2 + g4«4V 4 , 

h 

h = -00V2 (^iEi + ^ 2 E 2 ) x z 

e 

--p s V(p + ginivi + q 3 n 3 v 3 , 
n 



(G) 



where ctq = e 2 /2irh 

Tuning the magnetic field one can always fulfill the con- 
dition vt = 1 in average. But due to structural de- 
fects the local vt can be larger than unity in some ar- 
eas and smaller than unity in other areas. For simplic- 
ity, we consider the system with two statistically dis- 
tributed areas of equal fraction with local vti > 1 and 
vti < 1 {{vti + vt 2 )/2 = 1). In these areas the vor- 
tex densities for certain k are nonzero: n\ = n 2 — n v 
in vt < 1 areas, and n 3 = = n v in vt > 1 areas 
(n v = n {vTi - 1) = n Q (l - VT2))- 

We start from the analysis of transport properties of bal- 
anced bilayers (v\ = v 2 ). In this case the electrical charges 
of the vortices are equal in modulus. It is reasonable to 
assume that the viscosity coefficients are the same for all 
species of vortices (771 = n 2 = r/ 3 = 774 = n) . 

Solving Eqs.Q and (j4]) and substituting the solition into 
Eq. ([6]) we obtain 
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In the upper(lower) sign corresponds to vt > 1 
(i*r < 1) areas. Here we introduce the notation j3 — ithno, 
j± = ji ± j 2 and E± = Ei ± E 2 . The system considered 
is a two-component isotropic conducting medium. Each 
component is characterized by the parameters er+ = cr xxl 
a a = a ty (the parallel current conductivities) and o~ = 
a xx> a a — a xy (the counterflow conductivities). 

The exact expressions for the effective conductivities can 
be obtained by the method developed in [2^, [26| . In case of 
equal fractions of the components the effective quantities 
read as 
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(8) 
c s i + a s2 

in the expression for (cr^o) 



s(a)n ' n the expression for (c"^), the index n = 
1, 2 numerates the components. 
Using Eqs. Q and ([8]), we obtain 
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-, (a- ) = 0. (10) 



The effective parallel current and counterflow resistiv- 
ities are obtained from ([§)), (TTU|) by the operation of in- 
version of the conductivity tensor. We take into account 
that n v /no is the small parameter. In the leading order in 
n v /no the resistivities equal to 
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The negative sign of p x r x a9 means that the induced drag 
voltage is opposite to the voltage drop in the drive layer. 
Note that our model predicts rather small p x r x a9 - it is 
quadratic in the small parameter n v /nQ. 

One can see that for large viscosity 77 ^> (3 the resistiv- 



ities plx « p£f 



p x r x ve tx exp(— A/T), in an qualitative 



agreement with the experiment. The Hall resistivities p xy , 
Pxy"" 9 a nd Pxy lve are in the quantitative agreement with 
experimental data. 

Let us now consider the case of imbalanced bilayers. In 
imbalanced systems the vortex charges differ not only in 
sign but in absolute value. The viscosity is caused by the 
interaction of vortices with the pinning centers, so one can 
assume that the corresponding activation energy depends 
on the vortex charge. According to Table [TJ we introduce 
two viscocity parameters n a = r\\ — 774 and r\ = rj 2 — 773. 
For simplicity, we restrict the analysis to the case of large 
viscosities ^a^b 3> P- In this case one can neglect the 
difference of in vt > 1 and vt < 1 areas. Then, the 
relations between the currents and the fields are found to 
be 
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The inverse relations in linear in n v / no order are 
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It follows from (fT3"| that 
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(15) 
(16) 



Here p n , a $ notifies the resistivity in n-th layer. 

One can see that the longitudinal resistivities of differ- 
ent layers are determined by different viscosities. There- 
fore, the resistivities may demonstrate different tempera- 
ture dependences. It was observed experimentally that the 
activation energy for the longitudinal resistivities is higher 
in the layer with larger electron concentration. It corre- 
sponds to the grows of A K under increase of the vortex 
charge q K . 

In addition, we note that in the leading order the Hall 
resistivities (IT5|) are not sensitive to the imbalance as it 
was seen in experiments. 



defect areas work as a sort of an exclude volume. Thus, 
the effective conductivity a xx decreases and the effective 
resistivity p xx increases. On the other hand, if the con- 
ductivity a xx 9le is of the same order as the conductivity 
determined by the direct contribution of charged vortices 
the defect areas have only a small impact on the resistivity 



To describe the influence of defect areas quantitati vely 
we use the approach developed in 



26|. According to [26 1 



in general case of unequal fractions of the components 
the effective longitudinal and Hall conductivities of a two- 
component systems are given by the expressions 
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4. Possible origin of longitudinal drag resistivity 



D = A(l - />, A]Ki + (f 2 [p, X] - X 2 )a s 



We have shown that due to specific relation between 
the vortex charge and the filling factor v the resistivities 
p\ xx and Pi F x are very close to each other. They almost 
cancelled each other in the expression for the longitudinal 
drag resistivity, and the latter quantity (see Eq. ([TT])) is 
very small, in difference with experimental data that show 
considerable negative resistivity. 

It is important to note that in case of perfect exciton su- 
perfluidity the longitudinal drag resistivity should be pos- 
itive and equal to the drive resistivity. Indeed, any normal 
conductivity channel contributes into a xx , so the resistiv- 
ity plx should be nonzero. In the counterflow geometry the 
normal channel is shunted by superfluid excitons and can- 
not result in nonzero p£f . It yields pf x ag = p d J x lve = pL/2 

In case of slightly imperfect exciton superfluidity the 
situation might be the following. An additional normal 
conductivity channel (not connected with free vortices) in- 
creases both g\x and <j xx . It results in an increase of p\ x , 
but in a decrease p xx , so, it gives positive p x r x ag '. In partic- 
ular, bound vortex pairs (which are charged are carry zero 
vorticity) should work in this direction. Thus, the ques- 
tion on the origin of negative longitudinal drag remains 
open. 

In this section we will show that the negative drag can 
be accounted for the presence of defect areas, where the 
intcrlayer phase coherence is suppressed completely. For 
simplicity, we consider the case of balanced bilayers. 

On a qualitative level the impact of defect areas can 
be understood as follows. In these areas the counterflow 
and parallel current conductivities are equal each other 
and coincide (at least, approximately) with the conductiv- 
ities of a single layer with the filling factor v = 1/2. The 
longitudinal conductivity of a single layer <J xx gle is much 
smaller than the counterflow conductivity determined by 
the exciton channel, and in the counterflow geometry the 



A = 



1 



4cr s icr s 2 



\J {a si + <7 s2 ) 2 + (p a \ - Va2) 2 

1 



-\JIp7l - CT s2 ) 2 + (<Jal - Val) 2 



p is the fraction of the component 1, and 
f\p,X] = -(\- P ) (1-A) 
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Eq. (|17p was obtained by mapping of the system with 
nonzero a xy to the system with zero Hall conductivity. 
The function f[p,x] approximates the effective conduc- 
tivity (a) of an isotropic two-dimensional two-component 
system with zero Hall conductivity: (a) — crif\p, o^/oi], 
where Ui is the longitudinal conductivity of the i-th com- 
ponent, and p is the fraction of the component 1. The 
explicit expression (I19|) w as obtained in the random resis- 
tor network approach [27J. 

The three-component system can be reduced to the two- 
component one in the following way. The defect areas are 
considered as the component 1. The rest areas with imper- 
fect exciton superfluidity are considered as the component 
2 with the conductivity equal to the effective conductivity 
of the system considered in the previous section. 

Expressing the conductivities in <7o units, we have 
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Figure 1: Longitudinal (a) and Hall resistivities (b) in units of <Jq 
versus the concentration of the areas with destroyed interlayer phase 
coherence. Solid curves - counterflow resistivities, dashed curves - 
drag resistivities. The longitudinal drag resistivity is given with re- 
versed sign. The parameters ot\ = 0.05 and «2 = 0.025 are used. 

The model contains two conductivity parameters a\ 1 ai 
(normally, ai < a\ <C 1) and the parameter p (the defect 
area fraction). 

The resistivity tensor is obtained by inversion of the 
effective conductivity tensor. The effective resistivities de- 
pend on p. Typical dependences are shown in Fig. 1. 

One can see that the longitudinal drag resistivity is 
nonzero (and negative). Other feature seen from Fig. 1 are 
the essential increase of the longitudinal and the appear- 
ance of Hall resistivity in the counterflow channel. The 
Hall drag resistivity is equal approximately to a^ 1 (that 
correspond to the experiment) if p less and not to close to 
1/2. ' 

5. Conclusion 

Free vortices are dangerous for the superfluidity. But 
even in two-dimensional superfluid systems nucleation of 
vortices does not exclude perfect superfluidity. Below the 
BKT transition the vortices of opposite vorticities bind in 
pairs. Motion of such pairs, in difference with motion of 
single vortices, does not result in phase slips. The question 
is why it is not the case for quantum Hall bilayers. 

In quantum Hall bilayers each elementary charge exci- 
tation produces two fractionally charged vortices that re- 
pulse due to Coulomb forces. There exists a critical charge 
excitation concentration above which a gas of unpaired 
vortices emerges. Due to renormalization of the superfluid 
stiffness the critical concentration can be rather small. If 
the condition vt — 1 is fulfilled in average and local total 
filling factor deviates from unity, the excess electrons or 



holes induce vortex excitations that suppress the superflu- 
idity. If one trusts in the explanation presented one could 
expect that concrete values of counterflow, parallel current 
and drag resistivities are not universal quantities and may 
vary from sample to sample. 

The problem of negative longitudinal drag resistivity 
forces to put forward the idea on the existence of areas 
with completely destructed interlayer coherence. This idea 
looks quite reasonable taking into account that real geom- 
etry of arms should result in essential edge effects (the 
presence of edge areas with destructed interlayer coher- 
ence). The same effect can be caused by structure defects 
and impurities in GaAs heterostuctures. 

Nevertheless, one can hope that genuine exciton super- 
fluidity can be realized in bilayer quantum Hall systems. 
The only problem that special requirements on purity and 
perfectness of the samples should be fulfilled. 
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